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Abstract 

The density matrix renormalization group ("DMRG") discovered by White 
has shown to be a powerful method to understand the properties of many one 
dimensional quantum systems. In the case where renormalization eventually 
converges to a fixed point we show that quantum states in the thermody- 
namic limit with periodic boundary conditions can be simply represented by 
a special type of product ground state with a natural description of Bloch 
states of elementary excitations that are spin-1 solitons. We then observe 
that these states can be rederived through a simple variational ansatz making 
no reference to a renormalization construction. The method is tested on the 
spin-1 Heisenberg model. 
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In the past two years the density matrix renormalization group ("DMRG") method 
has been extensively tested on 1-d spin systems and spectacular numerical accuracy of 
both ground state energies and elementary excitations have been obtained with modest 
numerical effort. [1-4] In contrast to many of these efforts we explore the nature of the 
DMRG construction using a relatively few number of basis states to keep the numerical 
calculations simple. We have then been able to describe the thermodynamic limit of the 
ground state and single particle excitations in a way that generalizes very simply to arbitrary 
numbers of states. Our work shows that the fixed point limit of the DMRG leads to an ansatz 
form for the ground state and elementary excited states which can be explored variationally 
and is fundamentally independent of a renormalization scheme. 

For definiteness, we focus on the spin-1 antiferromagnetic Heisenberg chain with bi- 
quadratic interactions 

JV-l 

H = ^ S n ■ S n+ i — fi{S n ■ S n+ i) 2 . 

n=0 

In the DMRG scheme, we work recursively with blocks representing say n sites numbered 
from the left on which a set of states {|a)} are defined. In principle there are 3 n quantum 
states in the block but we keep only a set of n s "important" states in our basis which we 
expect accurately describe most operators in the ground state. 

In the DMRG recursion the next spin to the right of the block is absorbed into a new 
block which now has the 3n s states of the product representation {|a) n _i(x)|s n )}. The DMRG 
now provides a method to construct the projection operator A n a, ^ ,Sn ^ which projects these 
states to a set of new basis states now representing the important n s ' states of the larger 
block. This is written 

\ot)n = Y, A n a ' {P '" n) \ S n)®\P)n-l W 

/3,S„ 

where s n is the z component of spin at site n. Steven White's DMRG algorithm is a 
particularly effective way to compute a good projection operator. 

Two crucial points follow. First we perform a trivial change in notation: A^'^[s n ] = 
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A n a, ^ ,Sn \ Second, we observe that if a recursion fixed point exists for A i.e. A n [s] — > A[s] 
as n — > oo, we find \a) n = 

J2 {A[s n ]A[s n _ 1 }...A[s 1 }) a ' fj \s n s n _ 1 ...s 1 ) <g> |/3) 

s n t s n — l • • • 

where \/3) represents some state far away from n. The form of a wave function homogeneous 
in the bulk of the chain is now clear: For every n s X n s matrix Q } we define 

\Q) = E H QA[s n ]A[s n _ 1 }...A[s 1 ] (2) 

Thus \Q) represents a state that is uniform in the bulk but a linear combination of boundary 
conditions defined by \a) n on the left and \/3) on the right. [5,6] The special case of Q = 1 
the identity matrix leads to a translationally independent state with periodic boundary 
conditions. Since this state turns out to be normalized we write |f) = f). For the "AKLT" 
f3 = —1/3 model our ground state ansatz is exact [7] as are the "matrix product" states. 
[8,9] 

Several important properties of A follow. First the projection should preserve orthonor- 
mal bases: 8 a ^ a i = (a'\a). Using the recursion formula Eq. 1 and the orthogonality of the 
local spin states and previous block states, we find 

8a,a< = E A* a '^'[s']A^[s](s'\ S )(P'\P) 

PP'ss' (3) 

= Z.(A[s]Ai[s])»>»' . 

Hence in matrix form we find ^ s A[s]A^[s] = 1. 

We would also like our trial ground state with periodic boundary conditions to be an 
eigenstate of parity V where 

V\1) = J2 H A[s n }...A[ Sl ] )V\s n ... Sl ) (4) 
= £>( A[s n ]...A[ Sl ] )\ Sl ...s n ). (5) 

A sufficient condition is that there exists an invertible matrix Q-p where 

Q v A[s] = sign[V](A[s]) T Q v (6) 



where A T denotes transpose. Then V\l) = (sign[P]) n \l) . This assertion is shown by 
inserting QpQp 1 into the trace in Eq. (5) , then commuting Qp through each matrix inside 
the trace and ultimately canceling the Q^, 1 again. 

To find Qp multiply each side of Eq. (6) by A^[s] and sum over s. Using Eq. (3) we find 
that 

or/ = sign[p] j2( AT i s ]p T Q T /(A T is]r p 

STV 

= sign[P] h2(A T [3]r T A^[s]\ Q T / 
= sign[P] J2(A T [s] <8> A[s]t^^Q T / 

STV 

Thus, Q-p, if it exists, is the eigenvalue of the operator 2 s (A T [s] <S> A[s]) with value ±1. In 
cases that we have looked at we have had no trouble finding Q-p. 

Our basis states should form a representation of total spin. Since adding an integral spin 
does not mix half odd or integral reps, the block can consist of reps that are all half-odd or all 
integer spin. Keeping approximately twelve states we have found that the representations of 
half integer spin gives far better numerical results, and in this case twelve states consisting of 
two spin 1/2 and two spin 3/2 in the total spin representations optimizes numerical accuracy. 

In the diagram in Fig. 1 we have labeled on the left the "old" representations uniquely 
by the ordinal number 7 with the total spin j. New reps 7' and total spin j 1 are shown 
similarly on the right. Implicit in the labeling of the old states is total j z which is not 
shown. There are thus four irreducible representations of total spin with a total of twelve 
"old" states \"f,j z ). 

When adding the new spin there are 36 "intermediate" states in the product represen- 
tation falling into 10 different irreducible representations. We project from these reps back 
down to 4 irreducible "new" reps 7'. The nonzero projection terms P 1 r< are indicated 
by lines in Fig. 1. There are sixteen nonzero projection terms which are related by the 
requirement that our new states are orthonormal. 
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FIG. 1. The construction of the block states is shown when 12 states are kept in the basis. Old 
reps are on the left and new reps on the right. Each line represents a nonzero projection P r ' r< of 
basis reps. 

Expressing this relation mathematically, we let 7 uniquely label a representation of total 
spin and .7(7) denote the value of total spin. Each old state is thus labeled by I7, ra), where 
m is the z-component of total spin. The new states are then given by 

\i', m ') = Y, pl ' ,1 h^(Y), m ') (7) 

where I7, j(7'), m') denotes the intermediate states formed by I7) ® \s) written using total- 
spin basis. Since we demand that the projection operators preserve total j and ra, these 
states can be explicitly constructed using the Clebsch-Gordan coefficients (ji, mi|(j, m)(l, s)) 
in the following form 

l7>i(7')> m ') = Z)0'(7')> m 'IO'(7)> m )( 1 > s )) (I 5 ) ® l7, m )) • 
Inserting this into Eq. (7) we find that 
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\i,m')= Y, A (7 '' m,) ' (7 ' m) W(|5) ® | 7 ,m)) 

s,(7,m) 

where 

A (VX),(7,m) [a] = i^O'CVWIOM™) 

Thus, although the projection matrices A contain 3 X 12 X 12 numbers there are in fact only 
the few degrees of freedom available in P 1 n . 

There are sixteen parameters in P 1 n . However demanding orthonormality of the basis 
states turn out to yield six constraints. Furthermore we fix the freedom of mixing arbitrarily 
the two spin 1/2 and likewise the spin 3/2, giving a total of only eight free parameters. [10] 
We can then use a variational principle for the energy to determine these parameters. At 
this point it is clear that the DMRG plays no essential role in the construction aside from 
providing a guide which reps to keep. 

Expectation values of operators such as the Hamiltonian or correlations are given by 
sums of terms of the form 

UH1) = £{ SJ }{ S /}M A*[s n ']...A*[ Sl '] ) x 

tr{ A[s n }...A[ Sl } )(s n '...s 1 '\h\s n ...s 1 ). (8) 

To simplify this, we denote the identity for tensor products: (A ® B)^ a, ^^ r,v ^ = A a,T B^ i,v 
and the relation tr( B ) tr( A ) = tr{ B ® A ). We define the mapping M from 3x3 spin 
matrices M to n 2 s X n 2 s matrices by 

M = J2M s >, s A*[s'] ® A[s] 

s' s 

We let S = (S X} S y , S z ) define the spin-1 representation of total spin, and 1 = Id the 
"hat" mapping of the 3x3 identity matrix. The operator e t7TlSz is identified as the string 
operator. [11] 

Using Eq. (8), we see that expectation values of spin, energy and spin correlation function 
are given respectively by 
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(ll^i) = tr( r^s ) 

(11^- • S J+1 \1) = tr( i n ~ 2 ss ) 

(ll^- • s J+m \i) = tr( r-^sl^s ). 

while the string order parameter correlation function is given by 

(l|S'oe^ 5l k..e^ 5 ™- 1 )*S' m |l) = 

We note that these formulas are identical to those for one dimensional finite temperature 
classical statistical mechanical models, where the matrix 1 is identified as the transfer matrix 
and Sj the ordinary spin operator. 

The spectrum of correlation lengths, i.e. the collection of all possible exponential decay 
lengths £ of all correlation functions of the form (Oi(x)02(y)) — > ae~^ x ~ y ^^ is given by the 
negative of the logarithms of eigenvalues of 1 with a similar relation of e tvSz to the spectrum 
of string correlation lengths. We note that 1 is guaranteed to have an eigenvalue of 1 due to 
Eq. (3). Since the spin operator is orthogonal to this eigenvalue the next leading eigenvalue 
determines the decay of spin correlations. The eigenvalue of 1 in e tvSz gives the long range 
order in string correlations. [12] 

An ansatz for a Bloch state \Q,k) of momentum k is given by \Q, k) n = 

J2e mk tr( A[s n ]...A[s m+1 }QA[s m }...A[s 1 ] )\s n ... Sl ) 

m 

where Q is an arbitrary n s X n s matrix. Using the cyclicity of the trace, we find that 

n-l 

(£', k\Q, k) n = J2 H (Q ® i)(e T T(i ® Q')r- m ) (9) 

m=0 

The Hamiltonian HQ/ t Q(k,n) = (Q 1 , k\H\Q , k) n can be calculated similarly. When n is a 
power of 2 these sums can be computed recursively with an effort of log(n). 

We define the discrete Laplace transform of a series {a n }^L by F(\) = J2^=o a n e ~ nX - 
The discrete Laplace transform of Eq. (9) and of the Hamiltonian as a function of n is given 
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by a convolution product. Given the numerical values of 1 and e tvSz we can then obtain an 
analytic expression for the discrete Laplace transform and extract the leading behavior of 
a n as n — >■ oo. Defining HQ l)Q (k,n) = (Q' , k\H\Q , k) n and GQ' t Q(k,n) = (Q',k\Q,k) n it can 
simply be shown that when k ^ for large n 

H(k, n) = nH^k) + H (k) + 0(z) n 
<3(&,re) = Go(k) + O(z) n 

where \z\ ~ .8 < 1 represents the next-leading eigenvalue of the matrix 1. We can then solve 
for the ground state and excitation spectrum for a fixed value of k in the thermodynamic 
limit if the following equation can be solved: 

(nH 1 (k) + H (k))\Q, k) n = (nE + A k )G (k)\Q } k) n . 

for all n in the limit of large n. It can be shown that H\(k) oc Go(k) } which allows us to 
identify the ground state energy E per site through the proportionality H\(k) = E Go(k). 
The excitation spectrum is then given by 

H (k)\Q,k) n = A k G (k)\Q,k) n 

Similar formulas can be obtained for k = 0. We note that the last formula is a Hamiltonian 
for the excitation spectrum which makes no explicit reference to a ground state. 

These calculations have all been tested on the spin-1 Heisenberg chain, keeping the twelve 
states discussed before. The resultant 8-parameter family of trial ground states was explored 
to find the state of lowest energy for f3 in the range to 1. For f3 = the variational 
ground state energy was found to be E (f3 = 0) = —1.40138 , E (f3 = .6) = —2.9184 
and E (f3 = 1) = —3.98455. The most accurate ground state calculations indicate that 
E (P = 0) = -1.401484038971(4) [1,3] 

The result for f3 = 1 is to be compared with the exact value of —4 obtained in Ref. [13]. 
The parity operator has been computed in all cases and it is verified that the ground state 
has parity — 1^ where N is the number of sites. 
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FIG. 2. The spectrum for (3 = is shown. The lowest single particle triplet is shown as a solid 
line, with the lightly shaded region representing two-particle excitations and the dark region three 
particle excitations. Solid lines define the boundaries to the two and three particle continuum. 
Dotted lines indicate the spectrum of higher energy single-soliton excitations. The spin of these 
dotted excitations are, in order of increasing energy at k = it: 0,1,2,2,3,1,1,0. 

An important issue is whether or not we have a good ansatz form for the excitations. For 
the pure Heisenberg chain f3 = 0, we find the single particle spectrum shown in Fig. 2. The 
low lying triplet branch defines the gap A, = .4094, which is very good compared to the 
most accurately known [1,3,14,15] result of 0.410502(1) . Furthermore, we compute the spin 
wave velocity v = 2.452 to be compared to the calculations in Ref. [3] where v = 2.49(1) was 
obtained. Clearly we reproduce the single particle triplet excitations with high accuracy 
considering the few number of states in our basis. Our calculation also yields a detailed 
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spectrum of the lowest lying "single soliton" excitations shown by dotted lines in Fig. 2. 

Our second lowest energy excitation at k = it is a singlet shown by a dotted line in Fig. 2 
with A v (singlet) = 2.348. As a function of k, the second lowest single-particle excitation is 
either a singlet or a spin-2 object, as has also been observed in exact finite size calculations. 
[16] For the string order parameter, we find g(oo) = —.3759 whereas best estimates are 
[1] g(oo) = —.374325096(2). We find an asymptotic spin correlation length of / = 3.963 
compared to best estimates of / = 6.03. The severe truncation of our basis to only twelve 
states has resulted in the asymptotic correlations being quite poor, although we have verified 
that intermediate length spin-spin correlations are consistent with more precise calculations. 
[17] Parity of each of the elementary excitations is verified by checking the relation Eq. (6) 
with Q rather than the matrices A. 

The boundary to two particle excitations at a given value of k is shown in Fig. 2 computed 
explicitly by minimizing the sum of energies of excitations whose pseudomomentum sum to 
k, and similarly for the three particle excitations. These results are shown by the light 
shaded and dark shaded regions in Fig. 2. The picture fits well in with previously obtained 
results. 

We have similarly computed spectra for various values of f3. [18-20] Near f3 = .6 the 
excitation spectrum at k = it crosses zero and becomes negative! Our interpretation of 
this is that our ground state ansatz is deficient, and this shows up as a condensation of 
elementary excitations. It is to be noted that Oitmaa et al. [16] also found that numerically 
the gap appeared to vanish rapidly near this value of /3, although they too were unwilling 
to conclude that this persisted in the thermodynamic limit. 

Our calculations are consistent with two possible scenarios of what happens near f3 = .6. 
A special value of f3 could exist where the gap closes and signals a new phase. Or, the gap 
is in fact small and persists all the way to f3 = 1 but we do not see it due to our restricted 
ansatz for the ground state. In this case, we believe more accurate DMRG calculations will 
also have similar difficulties. 

A significant issue appears to be that the DMRG fixed point seems to invariably lead to a 
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matrix product ground state that, although it succeeds in reproducing ground state energies 
to high accuracy, cannot strictly give a power law decay of spin correlations. Thus, we find 
the energy very accurately at the Bethe ansatz point (3 = 1 without finding the expected 
powerlaw decay of correlations. The correlation length spectrum is given by the eigenvalues 
of the matrix 1, and it is hard to see how this can ever give algebraic correlations. However, 
intermediate correlations for intermediate lengths appear to be well represented in all cases. 
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